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ABSTRACT 

We report on a systematic investigation of the H/3 and Fe II emission lines in a sample of 568 quasars 
within z < 0.8 selected from the Sloan Digital Sky Survey. The conventional broad H/3 emission line can be 
decomposed into two components — one with intermediate velocity width and another with very broad width. 
The velocity shift and equivalent width of the intermediate-width component do not correlate with those of the 
very broad component of H/3, but its velocity shift and width do resemble Fe II, Moreover, the width of the very 
broad component is roughly 2.5 times that of the intermediate-width component. These characteristics strongly 
suggest the existence of an intermediate-line region, whose kinematics seem to be dominated by infall, located 
at the outer portion of the broad-line region. 

Subject headings: galaxies: nuclei — (galaxies:) quasars: emission lines — (galaxies:) quasars: general — 
galaxies: Seyfert — line: profiles 
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1. INTRODUCTION 

The geometry and kinematics of the broad-line region 
(BLR) in active galactic nuclei (AGNs) have been studied 
for about three decades but the details are far from well 
understood. It is widely accepted that the BLR is strati- 
fied: high-ionization lines originate from small radii and low- 
ioniza tion lines arise further out (Collin-Souffr in & Lasotal 
1988). This stratification picture is supported by the results 
of reverberation mapping, which show that lines of differen t 
ionization have different lags (e.g., Peterson & Wandel 1999). 
The dependence of the s ystemic velocities of the emission 
lines on ion ization (e.g., Gaskell 1982; Sulentic et al. 2000a; 
iRichards et al.l2002tlSn'ang et alJ20 07) sugg ests that the BLR 
may originate from a wind and a disk (e.g., Leighly & Moore 
l2004t iLeighlvl 12004. and references therein). However, the 
profiles of the broad emission lines often contain multiple ve- 
locity components, suggesting that the structure of the BLR 
may be more complex than can be described by a simple strat- 
ification or wind + disk model. 

It has been known that the broad H/3 line profiles are gen- 
erally not well described by a single G aussian. Two Gaus- 
sians (e.g., Netzer & Trakhtenbrot 2007) or a Gauss-Hermite 
function (e.g., Salviander et al. 2007) are often used. Addi- 
tionally, the profiles show great diversity from object to ob- 
ject. Sources with narrower H/3 lines tend to have stronger 
line wings, while those with broader H/3 lin es are dominated 
by the line core fe.g.. ISulentic et all 12002). Some sources 
have asymmetric and shifted H/3 profiles, suggesting that 
the BLR has a structure more complex than a single viri- 
alized component. The H/3 profile of OQ 208, for exam- 
ple, has an additional redshifted H/3 component of interme- 

1 National Astronomical Observatories of China, Chinese Academy of 
Sciences, Beijing 100012, China; chenhu@bao.ac.cn 

2 Key Laboratory for Particle Astrophysics, Institute of High Energy 
Physics, Chinese Academy of Sciences, Beijing 100049, China. 

3 Graduate University of the Chinese Academy of Sciences, Beijing 
100049, China. 

4 Theoretical Physics Center for Science Facilities (TPCSF), Chinese 
Academy of Sciences, China. 

5 The Observatories of the Carnegie Institution of Washington, 813 Santa 
Barbara Street, Pasadena, CA 91 101, USA. 

6 Department of Physics and Institute of Theoretical Physics, Nanjing Nor- 
mal University, Nanjing 210097, China. 



diate width that clos ely resembles the kinematics of Fe II 
(Marzi anTet al.lll993l) . Additional evidence that the BLR con- 
tains two or more kinematically distinct components comes 
from differential variability between the line core and wing 
(e.g jFerland et al] [l990: Peterson et al. 2000). The existence 
of an intermediate-line region and a very broad-line region for 
H/3 emission was suggested by some previous studies (e.g., 
lBrothertonll996llSulentic et alj|2000bl) . 

In a recent spectral decomposition of a large sample of 
quasars selected from the Sloan Digital Sky Survey (SDSS), 
Hu et al. (2008, hereinafter Paper I) find that the majority of 
quasars show Fe II emission that is both redshifted and nar- 
rower than H/3. Moreover, the magnitude of the Fe II redshift 
correlates inversely with the Eddington ratio. These charac- 
teristics suggest that Fe II originates from an exterior portion 
of the BLR, whose dynamics may be dominated by infall. 
These findings offer fresh insights into the structure of the 
BLR. 

In light of the trends associated with Fe II emission, we 
expect that a portion of the H/3-emitting gas may be related 
to an inflowing component too. In this Letter, we systemati- 
cally study the H/3 profiles of SDSS quasars to try to answer 
two questions: do all quasars have an intermediate-width H/3 
component similar to OQ 208, and, if so, is this component 
also associate with the Fe II emission? We show that the con- 
ventional broad H/3 line actually consists of two kinematically 
linked components, one of which originates from the same re- 
gion that emits Fe II. 

2. SAMPLE AND DATA ANALYSIS 

Our sa mple is selected from the S DSS Fifth Data 
Release dA delman-McCarthv et alJ l2007h quasar catalog 
dSchneider et alJ 120071) . We choose objects with redshifts 
z < 0.8 to ensure that [O III] A5007 lies within the SDSS spec- 
tral coverage. We also require a signal-to-noise ratio (S/N) > 
10 in the restframe wavelength range 4430-5550 A (cover- 
ing H/3, [O III], and the most prominent features of optical 
Fe II emission) and that no more than 1/3 of the pixels are 
masked by the SDSS pipeline in this region. 7601 quasars 
satisfy these criteria. The spectral analysis, whose details are 
described in Paper I, involves fitting a continuum model in a 
set of windows devoid of strong emission lines that consists of 
(1) a single power law, (2) Balmer continuum supplemented 
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FIG. 1 . — (a) Variation of H/3 profile with line width. P denotes the prob- 
ability that the double-Gaussian model cannot improve the fit. A high value 
of P means that the H/3 line can be adequately described by a single Gaus- 
sian, whereas a low value of P indicates that two Gaussians are needed. The 
dashed line is 1 -P = 99.73%. The sources below this line need two Gaus- 
sians to fit their broad H/3 lines, (b) The critical EW ratio of H/3ic/H/3vbc for 
requiring the double-Gaussian model as a function of signal-to-noise ratio. 

with high-order Balmer emission lines, and (3) a pseudo- 
continuum due to blended Fe II emission. The full width at 
half maximum (FWHM) and shift of Fe II are measured from 
the continuum decomposition. After subtracting the contin- 
uum model, the H/3 line is decomposed into narrow (H/3nc) 
and broad (H/3bc) components. H/3nc is forced to have the 
same profile as [O III], a shift of up to 600 km s" 1 relative to 
[O III], and an intensity constrained to lie between 1/20 and 
1/3 of tha t of [O III]. H/3rc is modeled u sing a Gauss-Hermite 
function (Ivan der Marel & Franxll 19931) . whose best fit yields 
FWHM(H/3bc)- The rest frame is defined by the peak of the 
[O III] A5007 (see Paper I for more details). 

Although the Gauss-Hermite function provides a good and 
convenient mathematical description of the H/3bc profile, for 
the purpose of this Letter — investigating whether there is 
an intermediate-width H/3 component associated with Fe II 
emission — we use the following simpler method to fit H/?bc- 
First, we use only one Gaussian to model the broad H/3 line 
(hereinafter the single-Gaussian model). Then we fit broad 
H/3 with two Gaussians (hereinafter the double-Gaussian 
model), an intermediate-width component (H/3ic) and a very 
broad component (H/?vbc)- In both models, H/3nc is fitted 
in the same way as in Paper I. First, we compare the \ 2 of 
the double-Gaussian model and the Gauss-Hermite. There 
are 630 sources with \ 2 20% larger in the double-Gaussia n 
model than in the Gauss-Hermite model (Hao et al. 2005); 
this means that the double-Gaussian fit is still inadequate for 
these sources 7 . Then, for the sources that can be well fitted by 
the double-Gaussian model, we compare the reduced x 2 of 
the do uble-Gaussian and single-Gaussian model, and use the 
F-test (Lupton 1993, Chapter 12.1) to calculate how signifi- 
cantly the double-Gaussian model improves the fit. Whether a 
source needs two Gaussians or not depends on its H/3bc width. 
Figure QJ shows the probability P that the double-Gaussian 
model cannot improve the fit, as a function of FWHM(H/3bc)- 
The 2435 sources below the dashed line (~ 32% of the whole 
sample) require two Gaussians at a significance greater than 
3 a. For sources with FWHM(H/3 BC ) > 5000 km s _I , one 
Gaussian is eno ugh to descr ibe the profile of H/3bc- This is 
consistent with [Collin et"ail (120061) . who find that the ratio 
FWHM/er (er is the line dispersion) of H/3bc depends on its 

7 From a statistical view for the whole sample, however, this does not mean 
that the double-Gaussian model is worse than the Gauss-Hermite model. 
There are 510 sources with x 2 20% smaller for the double-Gaussian model 
than for the Gauss-Hermite model, and for most sources both models are 
equally good. 



line width (see their Fig. 3). For sources with broader H/3bo 
this ratio is close to 2.35, the value for a single Gaussian. As 
discussed later, this means that broader sources have a weaker 
H/3ic component. 

We use Monte Carlo simulations to determine the detec- 
tion threshold of H/3ic- We generate artificial H/3 lines using 
two components [FWHM(H/3 V bc) is set to 2.5 FWHM(H/3ic) 
as shown in Fig. |5), fit it using a single and a double Gaus- 
sian, and then calculate the probability P as before. The EW 
ratio of H/3ic/H/3vbc is increased progressively until \—P ex- 
ceeds 99.73%. (All EW measurements refer to the continuum 
at 5100 A.) The critical EW ratio of H/3i C /H/3 V bc depends 
on S/N, as shown in Figure Q}?. In order to detect an H/3ic 
component with an EW that is 10% larger than the EW of 
H/3vbc, the S/N should be larger than 22.5. Thus, in the anal- 
ysis below, we remove all sources with (1) S/N < 22.5, (2) \ 2 
of double-Gaussian 20% larger than the \ 2 of Gauss-Hermite 
and 1 -P < 99.73%, and (3) EW(Fe) < 25 A (to insure that 
the Fe II measurements are as reliable as those in Paper I). 
These three cuts left 1499, 811, and 568 sources, respectively. 
Note that the fraction of sources that require a double-Gussian 
model (81 1/1499 sw 54%) is larger than the fraction when the 
S/N threshold is 10 (2435/7601 « 32%), consistent with the 
simulations above: fainter H/3ic can be detected with higher 
S/N. The analysis below focuses on the most stringent subset 
of 568 sources (~ 38% of the sources with S/N > 22.5). The 
reduced \ 2 f° r the double-Gaussian model fit has a median 
value of 1.004. 

To constrain the relative strengths of Fe II emission associ- 
ated with the very broad and intermediate-width components, 
we generated simulated spectra using two Fe II components, 
and then fit them using only one component, whose width is 
fixed to that of the input intermediate-width component. The 
EW of the input very broad component is increased progres- 
sively until the reduced \ 2 of the fit exceeds 1 a from the 
expect value. For a typical Fe II EW of 75 A and S/N sw 25 in 
our final sample, the flux ratio of Fe II between the input very 
broad component and intermediate-width component is ~0.3. 
This is the upper limit of the Fe II emission coming from the 
very broad component. 

Figure |2] shows the emission-line fitting and Fe II emission 
measurement for three typical sources. Example a is an ex- 
treme case that has a large Fe II redshift (1533 ±24 km s _1 ) 
and H/3 with an isolated red peak. The Fe II in example b 
is moderately redshifted (590±98 km s" 1 ), and its H/3 core 
shows only a strong red asymmetry rather than another peak. 
A more common situation is seen in example c, in which Fe 
II and H/3ic have no shift and the H/3 profile is symmetrical. 
Table Q] lists the measurements of the three sources shown in 
Figure 

3. RESULTS AND DISCUSSIONS 

We find that H/3ic and Fe II emission have similar kinemat- 
ics. The similarity can be seen not only in individual sources, 
but also statistically for the whole sample. From examples a 
to c in Figure [2] the H/3 profiles change progressively while 
the Fe II shifts become lower and lower. Figure 02 shows a 
strong correlation between H/3ic and Fe II shifts. Pearson's 
correlation coefficient rp is 0.22, and the probability P of a 
chance correlation is < 1 x 10~ 5 . H/3ic and Fe II have ap- 
proximately the same shifts. The widths of H/3ic and Fe II 
are also well correlated and roughly equal (Fig. 0?); rp = 
0.48 and P < 1 x 10~ 5 . Except for some sources to the upper 
left of the solid lines whose errors are large, the majority of 
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FIG. 2.— Examples of Fe II measurement and emission-line fitting for (a) SDSS J094603.94+013923.6, (b) SDSS J092008.22+032245.4, and (c) SDSS 
J103859. 58+422742. 2. For each source, the top panel shows multi-Gaussian fitting for H/3 and [O III]. H/3 v bc is m green and Hftc is in magenta. The two blue 
dashed lines mark the rest-frame wavelength of H/3 and [O III] A5007. The magenta dotted line is the position of H/3 at the same velocity as Fe II. Note the 
consistency between the Hftc peak and the dotted line. H/3nc and [O III] are in blue. The red line is the sum of each component. The bottom panel shows the 
emission-line spectrum after subtracting the power-law continuum. The green portions of the spectrum denote the windows for fitting the Fe II emission, whose 
model is given in red. The blue dashed line marks the peak of Fe II A4924 at zero velocity shift. 
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(kms -1 ) 




094603.94+013923.6 


0.220 


60±1 


I543±55 


I533±24 


54±I 


1428±18 


I754±7 


125 ±1 


5730± 58 


400±25 


3.506 


092008.22+032245.4 


0.334 


45 ±2 


2548±251 


590±98 


27±2 


2281±103 


589±42 


42±2 


6924±332 


-195±109 


0.715 


103859.58+422742.2 


0.221 


72±1 


I297±40 


-10±18 


I3±l 


1259±7I 


-I6±26 


44±1 


3893±87 


-279±29 


1.361 



NOTE. — Table[T]is available in its entirety electronically. 
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FIG. 3. — Correlations between (a) Fe II and H/3ic shifts and (b) Fe II and 
H/?ic widths. The solid diagonal lines denote that Fe II and H/?ic have the 
same shifts and widths. The contours in panel (a) show the density of the 
data points. 
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Fig. 5.— Plots of (a) H/3 V bc vs. H/3 rc shifts and (b) H/3 V bc vs. H/3 1C 
widths. The two components have different shifts and widths. The solid line 
in panel (b) denotes FWHM(H/3 VB c) = 2.5 FWHM(H/3 IC ). 
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FIG. 4. — Plots of (a) Fe II vs. H/3 V bc shifts and (b) Fe II vs. H/3 V bc 
widths. The two lines have different shifts and widths. 



sources follow the relation that Hj3ic an d F e n nave me same 
shifts and widths. By contrast, Fe II and H/3vbc have different 
shifts and widths (Fig. |4j, We also find that the kinematic 
connection between H/3ic and H/3vbc is complicated. No cor- 
relation between H/3rc an d H/3vbc shifts is seen (Fig. [5}0, 
but the H/3ic an d H/3vbc widths are strongly linked, such that 
FWHM(H/3 V bc) « 2.5 FWHM(H/3ic) (Fig. |3J>). 

The observational results described above suggest a sce- 
nario in which the conventional BLR consists of two 
components — an intermediate-line region (ILR) and a very 
broad-line region (VBLR). If both regions are virialized, so 
that R oc v~ 2 , then the ILR is about 2.5 2 or 6.25 times farther 
from the center than the VBLR. Because of its redshift, the 
kinematics of the ILR may be dominated by infall. H/3 emis- 
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FIG. 6.— Plots of (a) H/3 V bc vs. H/3 IC EWs and (b) Fe II vs. H/3 IC EWs. 
No correlations are found. 

sion emerges from both the ILR and VBLR, while most of the 
Fe II emission comes from the ILR. 

The ILR and VBLR defined in this paper a re es - 
sentially s i milar to those descr i bed in ICorbinl (1 19951). 
Brotherto^ (119961) . ISulentic et afl (2000b), and IZhuetalJ 
(120081) but differ from tho se in |Brotherton ~et"aIl d 19941) or 
ISulentic & Marzianil (11999b . which refer to the C iv-emitting 
region. Note that the ILR and VBLR of the H/3-emitting re- 
gion are distinct from those of the C iv-emitting region be- 
cause the C IV ILR usually has the systemic redshift while 
the C IV VBLR shows a blueshift. 

The lack of correlation between the EWs of H/3ic and 
H/?vbc (Fig. |6}i) strongly suggests that the two components 
are emitted from different regions. If both are photoion- 
ized, they must have different covering factors. The relative 
strength between H/3ic and H/3vbc determines the final H/3 
profiles. In FigureQ] sources broader than ^5000 km s" 1 can 
be fitted well using one Gaussian. This trend can be inter- 
preted in our two-component BLR scenario. It is well known 
that sources with broad H/3 tend to have weak Fe n/H/3 (e.g., 
lBoroson&Greenlll992t ISulentic etai1l2000al) . The ILR is 
weak in these systems because Fe II is weak. Their profiles 
show little deviation from a single Gaussian under the typical 
S/N level of SDSS spectra. The composite spectra of sources 
with large Fe II redshifts, on the other hand, do show red 
asymmetry in the H/3 profiles (see Fig. 13 of Paper I). The 
variation in the relative strength of the H/3vbc and Hftic com- 
ponents in different sources reflects the competition between 
the two components, although they apparently do so in such a 
manner that their kinematics remained coupled. 

It is of interest to note that the strengths of H/3ic and Fe II 
are not correlated (Fig. |6j>). The wide range of Fe u/Hftic 



ratios reflects either the comp lexity of the e xcitation mecha- 
nism of Fe II emission (e.g., Bald win et al.l 12004. and refer- 
ence^ therein) or lar ge variat ions in quasar metallicities (e.g., 
iNetzer & Trakhtenbrotl2007h . 

4. SUMMARY 

We have studied the profiles of the H/3 emission line using 
a large sample of quasars selected from SDSS. Comparing the 
H/3 profiles with the properties of Fe II emission given in Pa- 
per I, we deduced the existence of two H/3 emission regions — 
an intermediate-line region and a very broad-line region. The 
observational evidence can be summarized as follows: 

1 . The velocity shifts and widths of the H/3 intermediate- 
width component are approximately the same as those 
of Fe II, indicating that they originate from the same re- 
gion. However, the Fe ll/H/3ic ratios vary greatly from 
object to object, reflecting variations in either metal 
abundance or excitation conditions in the Fe Il-emitting 
region. 

2. The velocity width of the very broad component of H/3 
is roughly 2.5 times larger than that of the intermediate- 
width component, but no correlation exists between 
their radial velocities. The equivalent widths of the 
two components are also unrelated, suggesting that 
they have very different covering factors and geome- 
try. The conventional BLR seems to consist of two dif- 
ferent, physically distinct regions. We suggest that the 
intermediate-width component of H/3 and Fe II trace an 
infalling region in the outskirts of the BLR, likely lo- 
cated in between the molecular torus and the accretion 
disk. 

The properties of the H/3-emitting region discussed here, in 
conjunction with those of Fe II emission summarized in Pa- 
per I, offer important new constraints on models of the broad 
emission-line regions. This will be the subject of a forthcom- 
ing paper. 
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